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METHOD FOR CONTROLLING GROWTH OF A SILICON CRYSTAL 
TO MINIMIZE GROWTH RATE AND DIAMETER DEVIATIONS 



CROSS REFERENCE TO RELATED APPLICATION 

This application is a continuation of commonly assigned application Serial No. 
09/495,563, filed February 1, 2000, the entire disclosure of which is incorporated herein by 
reference. 



BACKGROUND OF THE INVENTION 

This invention relates generally to improvements in controlling growth processes 
of single crystal semiconductors for use in the manufacture of electronic components and, 

10 particularly, to methods of accurately controlling growth in a Czochralski crystal growth 

process for minimizing deviations in both growth rate and diameter, including determining 
steady-state growth rate and controlling pull rate around the steady-state growth rate, 
independent of process conditions such as sensed melt temperature. 

Monocrystalline, or single crystal, silicon is the starting material in most processes 

15 for fabricating semiconductor electronic components. Crystal pulling machines employing 

the Czochralski process produce the majority of single crystal silicon. Briefly described, 
the Czochralski process involves melting a charge of high-purity polycrystalline silicon in 
a quartz crucible located in a specifically designed furnace. After the heated crucible melts 
the silicon charge, a crystal lifting mechanism lowers a seed crystal into contact with the 

20 molten silicon. The mechanism then withdraws the seed to pull a growing crystal from the 

silicon melt. A typical crystal lifting mechanism suspends the seed crystal from one end of 
a cable, the other end of which is wrapped around a drum. As the drum rotates, the seed 
crystal moves up or down depending on the direction that the drum is rotating. 

After formation of a crystal neck, the growth process enlarges the diameter of the 

25 growing crystal by decreasing the pulling rate and/or the melt temperature until a desired 

diameter is reached. By controlling the pull rate and the melt temperature while 
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compensating for the decreasing melt level, the main body of the crystal is grown so that it 
has an approximately constant diameter (i.e., it is generally cylindrical). Near the end of 
the growth process but before the crucible is emptied of molten silicon, the process 
gradually reduces the crystal diameter to form an end cone. Typically, the end cone is 
formed by increasing the crystal pull rate and heat supplied to the crucible. When the 
diameter becomes small enough, the crystal is then separated from the melt. During the 
growth process, the crucible rotates the melt in one direction and the crystal lifting 
mechanism rotates its pulling cable, or shaft, along with the seed and the crystal, in an 
opposite direction. 

Although presently available Czochralski growth processes have been satisfactory 
for growing single crystal silicon useful in a wide variety of applications, further 
improvements are still desired. For example, conventional growth processes employ one 
of two main strategies for crystal growth control. First, standard control aims primarily to 
control diameter of the growing crystal. Second, "locked seed lift" control aims primarily 
to control growth rate. It is desired to minimize deviations in both diameter and growth 
rate simultaneously throughout the growth process to improve the shape and the quality of 
the as-grown crystal. 

The problem of minimizing growth rate variations involves reliably measuring 
growth rate. In conventional growth processes, the pull rate, which is set by the seed lift, 
is a known parameter. The growth rate can be determined indirectly by examining the 
meniscus dynamics if the height of the meniscus at the crystal-melt interface is known as a 
function of time. In an ideal case, the ingot has a perfectly cylindrical shape and the pull 
rate is identical to the growth rate. Achieving relatively small deviations in crystal 
diameter is coupled to achieving relatively small deviations in growth rate from the 
prescribed value. Although attempts have been made using pyrometers or laser reflections 
from meniscus shape, a reliable system for measuring meniscus height is not available. 
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Therefore, improvements in crystal growing processes are desired for determining growth 
rate without measuring meniscus height. 

A controller is also desired for solving the problem of unstable crystal grow th. In 
particular, such a controller is needed for maintaining a stable crystal growth for the 
duration of the growth process by minimizing diameter variations. Those skilled in the art 
recognize that pulling a crystal from a melt is an intrinsically unstable process. Without 
control, the crystal will grow in a conical, rather than generally cylindrical, shape. Active 
control loops are necessary to maintain the cylindrical shape, which means keeping the 
pull rate close to the steady-state growth rate. Perturbations in the steady-state growth rate 
can be time-dependent and can differ from run to run and/or puller to puller. 

As described above, conventional growth processes use two main strategies for 
crystal growth control. The first primarily involves diameter control and the second 
primarily involves a kt locked seed lift" for growth rate control. The presently available 
systems, how ever, do not provide simultaneous control of diameter and growth rate. For 
example, standard diameter control manages crystal growth by means of two PID 
(Proportional Integral Derivative) control loops acting on pull rate and heater power. The 
primary objective of the controller is to grow a crystal at diameter equal to a set value 
(given other process parameters such as crystal and crucible rotation or the presence of 
magnetic fields). The first PID control unit acts on pull rate to minimize the diameter 
error. As a secondary objective, the puller controller seeks to minimize the error between 
the average pull rate and the set value. This is accomplished by means of the second PID 
control loop acting on heater power. With the standard PID control system, diameter 
variation can be set to a tolerance of about 1 mm but relatively large pull rate deviations 
from the set value can still occur (e.g., even of the order of 30-40%). 

The primary objective of conventional seed lift control systems is growth rate 
control. In this instance, pull rate is set to a constant prescribed value V sct . Diameter 
control is of secondary importance and is generally accomplished by means of a PID 
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control unit acting directly on the heater power. Although this control strategy precisely 
controls pull rate relative to the set value, relatively large diameter deviations (e.g., up to ± 
5-10 mm) are possible. Moreover, the actual growth rate still deviates even though the 
pull rate equals V sct so the main goal of growth rate control is not completely achieved. 

By properly tuning the controller, acceptable performance with respect to the 
growth rate may be possible with presently available control systems. In this instance, 
proper tuning means empirically setting the various controller coefficients by trial and 
error given a particular hot zone configuration and a set of process parameters. The set of 
coefficients is defined for different portions of the crystal length to compensate for 
changes in the thermal configuration caused by the batch nature of the process. Given the 
number of different zones taken into consideration and the type of controller, a very large 
number of coefficients (e.g., 50-100 or more) must be properly tuned to obtain acceptable 
growth characteristics. Any change in hot zone design or process conditions, such as 
crucible or crystal rotations or an applied magnetic field, can negatively impact the 
performance of the controller. More re-tuning iterations are then necessary to 
accommodate such changes. In addition, critical processes where minimum tolerances are 
set for diameter and growth rate deviations often must be personalized for each puller. 
Tuning and optimization of process control parameters, updating of control parameters to 
new hot zone configurations or process conditions, and management of puller dependent 
processes is very cumbersome and costly. Thus, a reduction in the number of iterations 
needed to achieve acceptable process performance or a reduction in resources dedicated to 
process maintenance and personalization will provide significant cost savings in producing 
superior silicon crystals. 

For these reasons, an accurate and reliable apparatus and method for controlling 
silicon crystal growth to minimize both grow th rate deviations and diameter deviations 
that is less costly and easier to tunc is desired. 
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SUMMARY OF THE INVENTION 

The invention meets the above needs and overcomes the deficiencies of the prior 
art by providing methods of minimizing both growth rate deviations and diameter 
deviations in a crystal ingot pulled from a melt according to the Czochralski process. 
Among the several objects of the invention may be noted the provision of such methods 
that reduce the number of control parameters; the provision of such methods that permit 
determining the control parameters independently of crystal position, hot zone 
configuration, process conditions and/or puller; the provision of such methods that permit 
determining the control parameters from observation of diameter evolution; the provision 
of such methods that predict desired interventions on the pull rate and the power; the 
provision of such methods that permit tuning the control parameters during the run; the 
provision of such methods that permit deriving the meniscus height dynamics from 
measured radius dynamics; the provision of such methods that permit determining 
meniscus height measurements as a function of the derived meniscus height dynamics; the 
provision of such methods that may be incorporated into existing crystal pulling devices; 
and the provision of such methods that can be carried out efficiently and economically. 

Briefly described, a method embodying aspects of the invention is for use with a 
crystal puller for growing a monocrystalline semiconductor crystal according to the 
Czochralski process. The crystal puller has a heated crucible containing a semiconductor 
melt from which the crystal is grown. The crystal puller also has a heater supplied by a 
power supply for heating the crucible. The crystal is grown on a seed crystal pulled from 
the melt at a pull rate V p . The method includes defining an initial interval of time for 
observing growth of the crystal being pulled from the melt and determining variations in 
crystal diameter occurring during the observation interval. Based on the variations in the 
crystal diameter, the method estimates a current growth rate V„ at which the crystal is 
being grown. The method also includes estimating a current steady-state value of the 
growth rate V, JS as a function of the estimated growth rate at the end of the observation 
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interval and determining a pull rate parameter and a heater power parameter as a function 
of the estimated steady-state growth rate V . The pull rate parameter is representative of 
an incremental change in the pull rate V to effect a desired change in the diameter of the 
crystal tow ard a target diameter and the heater power parameter is representative of an 
5 incremental change in the power supplied to the heater to effect a desired change in the 

growth rate of the crystal toward a target growth rate. According to the invention, the 
parameters are determined independently of a temperature condition sensed during pulling. 
The method further includes adjusting the pull rate V p according to the pull rate parameter 
and adjusting the power supplied to the heater by the power supply according to the heater 

10 power parameter. This minimizes variations in both crystal diameter and growth rate 

during subsequent growth of the crystal following the initial observation interval. 

Another embodiment of the invention is directed to a control method for use with a 
crystal puller for growing a monocrystalline semiconductor crystal according to the 
Czochralski process. The crystal puller has a heated crucible containing a semiconductor 

1 5 melt from which the crystal is grown. The crystal puller also has a heater supplied by a 

power supply for heating the crucible. The crystal is grown on a seed crystal pulled from 
the melt at a pull rate V p . The melt has a surface with a meniscus adjacent the crystal as 
the crystal is pulled from the melt. The method includes defining an initial interval of time 
for observing growth of the crystal being pulled from the melt and determining variations 

20 in crystal diameter occurring during the observation interval. In addition, the method 

includes defining a function r(t) based on the variations in crystal diameter occurring 
during the observation interval. The function r(t) represents radius variations and is a 
function of current values of crystal radius r, meniscus height h and growth rate V tt with 
respect to time. The method then performs a best fit routine on the function r(t) to deduce 

25 the current values of crystal radius r f , meniscus height h t and grow th rate V gf at the end of 

the observation interval. In turn, a current steady-state value of the growth rate V gs is 
calculated as a function of the deduced current values of crystal radius r t , meniscus height 
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h t and growth rate V f at the end of the observation interval and independent of melt 
temperature and meniscus height sensed during pulling. The method further includes 
controlling the crystal puller as function of current steady-state growth rate V gs to 
minimize variations in both crystal diameter and growth rate during subsequent growth of 
5 the crystal. 

In yet another embodiment, a method of the invention is for use with a crystal 
puller for growing a monocrystalline semiconductor crystal according to the Czochralski 
process. The crystal puller has a heated crucible containing a semiconductor melt from 
which the crystal is grown. The crystal puller also has a heater supplied by a power supply 

10 for heating the crucible. The crystal is grown on a seed crystal pulled from the melt at a 

pull rate V . The method includes defining an initial interval of time for observing growth 
of the crystal being pulled from the melt and determining variations in crystal diameter 
occurring during the observation interval. Based on the variations in the crystal diameter, 
the method estimates a current growth rate V g at which the crystal is being grown. The 

15 method also includes estimating a current steady-state value of the growth rate V gs as a 

function of the estimated growth rate V uf at the end of the observation interval and 
independent of meniscus height measured during pulling. In turn, the method determines a 
pull rate parameter and a heater pow er parameter as a function of the estimated steady- 
state growth rate V gs . The pull rate parameter is representative of an incremental change in 

20 the pull rate V p to effect a desired change in the diameter of the crystal toward a target 

diameter and the heater power parameter is representative of an incremental change in the 
power supplied to the heater to effect a desired change in the growth rate of the crystal 
toward a target growth rate. The method further includes adjusting the pull rate V p 
according to the pull rate parameter and adjusting the power supplied to the heater by the 

25 power supply according to the heater power parameter. This minimizes variations in both 

crystal diameter and growth rate during subsequent growth of the crystal following the 
initial observation interval. 
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Alternatively, the invention may comprise various other methods and systems. 
Other objects and features will be in part apparent and in part pointed out 
hereinafter. 

BRIEF DESCRIPTION OF THE DRAWINGS 
5 FIG. 1 is a partially schematic, block diagram illustrating a crystal puller and an 

apparatus according to the present invention for controlling the crystal puller. 

FIG. 2 is a schematic, fragmentary view of a silicon crystal being pulled from a 
melt contained in the crystal puller of FIG. 1. 

FIG. 3 is an exemplary flow diagram illustrating operation of the apparatus of FIG. 
10 1 to generate control signals for controlling the crystal puller. 

FIG. 4 is an exemplary graph illustrating crystal radius deviations as affected by 
incremental changes in pull rate and heater power as controlled by the apparatus of FIG. 1. 

FIG. 5 is an exemplary graph further illustrating crystal radius deviations as 
affected by incremental changes in pull rate and heater power as controlled by the 
15 apparatus of FIG. 1. 

Corresponding reference characters indicate corresponding parts throughout the 
drawings. 

DETAILED DESCRIPTION OF PREFERRED EMBODIMENTS 

Referring now to FIG. 1, a control apparatus, indicated generally at 1 1, is shown 

20 for use with a Czochralski crystal grow ing apparatus, indicated generally at 13. The 

details of construction of the crystal growing apparatus, or puller, 13 are well known to 
those of ordinary skill in the art. In general, the crystal puller 13 includes a vacuum 
chamber 15 enclosing a crucible 19. Heating means such as a resistance heater 21 
surrounds the crucible 19. In one embodiment, insulation 23 lines the inner wall of 

25 vacuum chamber 15 and a chamber cooling jacket (not shown) fed with water surrounds it. 
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A vacuum pump (not shown) typically removes gas from within the vacuum chamber 15 
as an inert atmosphere of argon gas is fed into it. 

According to the Czochralski single crystal growth process, a quantity of 
polycrystalline silicon, or polysilicon, is charged to crucible 19. A heater power supply 27 
5 provides electric current through the resistance heater 21 to melt the charge and, thus, form 

a silicon melt 29 from which a single crystal 31 is pulled. As is known in the art, the 
single crystal 3 1 starts with a seed crystal 35 attached to a pull shaft, or cable, 37. As 
shown in FIG. 1, single crystal 31 and crucible 19 generally have a common axis of 
symmetry 39. One end of cable 37 is connected by way of a pulley (not shown) to a drum 

10 (not shown) and the other end is connected to a chuck (not shown) that holds the seed 

crystal 35 and the crystal 31 grown from the seed crystal. Commonly assigned U.S. Patent 
No. 5,935,328, the entire disclosure of which is incorporated herein by reference, shows a 
drum and cable arrangement suitable for use with the present invention. 

During both heating and crystal pulling, a crucible drive unit 43 rotates crucible 19 

15 (e g-, in the clockwise direction). The crucible drive unit 43 also raises and lowers 

crucible 19 as desired during the growth process. For example, crucible drive unit 43 
raises crucible 19 as the melt 29 is depleted to maintain its level, indicated by reference 
character 45, at a desired height. A crystal drive unit 47 similarly rotates the cable 37 in a 
direction opposite the direction in which crucible drive unit 43 rotates crucible 19. In 

20 addition, the crystal drive unit 47 raises and lowers crystal 31 relative to the melt level 45 

as desired during the growth process. 

In one embodiment, crystal puller 13 preheats the seed crystal 35 by lowering it 
nearly into contact with the molten silicon of melt 29 contained by crucible 19. After 
preheating, crystal drive unit 47 continues to lower seed crystal 35 via cable 37 into 

25 contact with melt 29 at its melt level 45. As seed crystal 35 melts, crystal drive unit 47 

slowly withdraws, or pulls, it from the melt 29. Seed crystal 35 draws silicon from melt 
29 to produce a growth of silicon single crystal 31 as it is withdrawn. Crystal drive unit 47 
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rotates crystal 31 at a reference rate as it pulls crystal 31 from melt 29. Crucible drive unit 
43 similarly rotates crucible 19 at another reference rate, but usually in the opposite 
direction relative to crystal 3 1 . 

The details of construction of crystal drive unit 47 are well known to those of 
5 ordinary skill in the art. Generally, drive unit 47 includes a motor (not shown) with its 

shaft mechanically coupled to the drum. Although this mechanical coupling may comprise 
a direct connection between the shaft and drum, a preferred arrangement positions a set of 
reduction gears between the two for better control and smoother operation. The motor is 
thus operable to let out and reel in cable 37 via the drum for lowering seed crystal 35 into 

1 0 melt 29 and for pulling the ingot 3 1 from melt 29. 

A control unit 51 initially controls the withdrawal, or pull, rate and the power that 
pow er supply 27 provides to heater 21 to cause a neck down of crystal 3 1 . Preferably, 
crystal puller 13 grows the crystal neck at a substantially constant diameter as seed crystal 
35 is drawn from melt 29. For example, the control unit 51 maintains a substantially 

15 constant neck diameter of about five percent of the desired body diameter. Under a 

conventional control scheme, control unit 51 adjusts the rotation, pull and/or heating 
parameters after the neck reaches a desired length to cause the diameter of crystal 31 to 
increase in a cone-shaped manner until a desired crystal body diameter is reached. For 
example, control unit 51 decreases the pull rate to create an outwardly flaring region 

20 typically referred to as the taper of the crystal. Once the desired crystal diameter is 

reached, control unit 51 controls the growth parameters to maintain a relatively constant 
diameter as measured by system 1 1 until the process approaches its end. At that point, the 
pull rate and heating are usually increased for decreasing the diameter to form a tapered 
portion at the end of single crystal 31 . 

25 Commonly assigned U.S. Patent No. 5,178,720, the entire disclosure of which is 

incorporated herein by reference, discloses one preferred method for controlling crystal 
and crucible rotation rates as a function of the crystal diameter. Commonly assigned U.S. 
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Patent No. 5,882,402, U.S. Patent No. 5,846,318, U.S. Patent No. 5,665,159 and U.S. 
Patent No. 5,653,799, the entire disclosures of which are incorporated herein by reference, 
provide accurate and reliable measurements of a number of crystal growth parameters, 
including crystal diameter. In these patents, an image processor processes images of the 
crystal-melt interface to determine the diameter. Commonly assigned application Serial 
No. 09/372,897, filed August 12, 1999, the entire disclosure of which is incorporated 
herein by reference, teaches accurately pulling a single crystal silicon ingot from a melt 
according to a predetermined velocity profile, or target, specified in a crystal "recipe" for 
controlling the formation of defects. 

FIG. 2 illustrates an intermediate phase of the crystal growth process following 
melt-down and dipping of seed crystal 35. After formation of a crystal neck, the typical 
process enlarges the diameter of the growing crystal 31 by decreasing the pulling rate 
and/or the melt temperature until a desired diameter is reached. This portion of increasing 
diameter is referred to as a taper, or crown. As the taper increases to a desired diameter, 
crystal 31 forms a shoulder followed by a main body. As melt 29 approaches depletion, 
the crystal diameter is gradually decreased resulting in the formation of a generally conical 
end cone of crystal 31. Once the diameter of the end cone is sufficiently small (e.g., 2 mm 
to 4 mm), detachment of crystal 31 from melt 29 can be accomplished without causing 
dislocations to spread to the main body of crystal 31. Crystal 31 is then removed from 
vacuum chamber 15 for processing into wafers. By controlling the pull rate and the melt 
temperature while compensating for the decreasing melt level, the body of crystal 31 is 
grown so that it has an approximately constant diameter and constitutes a generally 
cylindrical body of crystalline silicon (i.e., an ingot). 

As the crystal drive unit 47 pulls crystal 31 from melt 29, a liquid meniscus 55 
forms at the interface between crystal 3 1 and melt 29. The liquid meniscus 55 forms on a 
surface 57 of melt 29. As is known in the art, the reflection of crucible 19 on meniscus 55 
is often visible as a bright ring adjacent crystal 31 . FIG. 2 also indicates the height h and 
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slope (J) of meniscus 55. In addition, FIG. 2 indicates the relative directions of the pulling 
velocity, or pull rate, V p and the growing velocity, or growth rate, V„. 

It should be understood that the as-grown crystal body typically will not have an 
entirely uniform diameter, although it is generally cylindrical. For this reason, the 
5 diameter of crystal 3 1 may vary slightly at different axial positions along axis 39. Further, 

the diameter varies in the different phases of crystal growth (i.e., seed, neck, taper, body 
and end cone). Crystal diameter measurements may be determined by a number of 
methods. For example, several technologies are known for providing crystal diameter 
measurements including methods of measuring the width of the bright ring. The bright 

10 ring is a characteristic of the reflection of the crucible wall in the meniscus that is formed 

at the solid-liquid interface of crystal 3 1 . Conventional bright ring and meniscus sensors 
employ optical pyrometers, photocells, rotating mirrors with photocells, light sources with 
photocells, line-scan cameras, and two-dimensional array cameras. Commonly assigned 
U.S. Patent Nos. 5,665,159 and 5,653,799, the entire disclosures of which are incorporated 

1 5 herein by reference, describe a system and method, respectively, for accurately and reliably 

measuring crystal diameter for use in controlling the grow th process of single crystal 
silicon. Advantageously, the system and method of these patents accurately determine the 
grow ing crystal's diameter by processing images of the crystal-melt interface generated by 
a camera. 

20 Referring again to FIG. 1, control unit 51 includes a programmable logic controller 

(PLC) 61 having a central processing unit (CPU) 63 and a memory 65. In one 
embodiment, the PLC 61 receives a position signal representative of the rotational position 
of the drum and is readily programmed via conventional means to convert the position 
signal into a numeric figure representing the linear movement of cable 37 in real time. In 

25 this manner, PLC 61 calculates the relative position of seed crystal 35 and the pull velocity 

of cable 37 as a function of the position signal. 
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In a preferred embodiment of the invention, the PLC 61 of control unit 51 
constitutes a model-based controller designed to achieve proper control of the growth rate 
and diameter of the growing crystal 31. Control unit 51 determines the current growth rate 
V u , as well as the current steady-state growth rate V gs , from radius dynamics. Control unit 
5 51 also operates to control the pull rate around V us to minimize the deviation in diameter 

from a set value. As described below, control unit 51 determines a preferred set of control 
parameters for the next fraction of growth by analyzing diameter excursions previously 
measured in a given length of grown crystal 31 . Such a diameter profile for a given length 
of crystal 31 provides information regarding radius variations occurring during an 

10 observation interval. 

The present invention advantageously controls silicon crystal growth to minimize 
growth rate and diameter variations independent of specific process conditions such as 
sensed melt temperature. The problem of minimizing growth rate variations involves 
reliably measuring growth rate. The invention determines a preferred set of control 

15 parameters for the next fraction of growth, following the observation interval, by analyzing 

diameter excursions previously measured in a given length of grown crystal. Such a 
diameter profile for a given length of crystal provides information regarding radius 
variations occurring during the observation interval, which is in contrast to merely taking 
diameter measurements at set intervals (i.e., periodically measuring or sampling the 

20 diameter). 

The control unit 51 preferably uses a reduced set of control parameters including 
two basic response constants that are a function of crystal length (e.g., a height response 
coefficient A h and a power response coefficient A p ). A third parameter, namely, transient 
time t, is less important in this embodiment. Control unit 51 preferably deduces the 

25 parameters A h and A from a single control run using the series of power pulses or from a 

global simulation of the temperature field in puller 13 for several values of the meniscus 
height and power. 
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Minimizing growth rate variations involves first reliably measuring growth rate. In 
conventional growth processes, the pull rate V p is a known parameter since it is set by the 
seed lift recipe. This recipe also includes a unique height (denoted by h s ) for growing a 
cylindrical ingot of a fixed diameter and any deviation of the meniscus height h from h s 
5 results in a diameter change. In other words, an evolution in diameter occurs because of 

an error in the meniscus height as a function of time, i.e., h(t). Equation (1 ), below, 
defines the radius deviation as a function of the meniscus height deviation: 

ar(t)/ct = - (2h s V g (t)/a 2 ) (h(t) - h s ) ( 1 ) 

where h(t) is the meniscus height (mm); h s is the steady-state meniscus height {e.g., 
10 h s = 6.5 mm); r(t) is the crystal radius (mm); r s is the crystal radius set point (mm); V g is 

the growth rate; and a is the capillarity constant (e.g., a = 7.6 mm). 

The meniscus height evolves if the pull rate V p deviates from the actual growth rate 
V a . Thus, if the height h(t) of the meniscus at the interface between crystal 31 and melt 29 
is known, the growth rate can be determined indirectly by equation (2), defined below: 

15 dh(t)/di = V p (t) - V g (t) (2) 

In an ideal case, the body of crystal 3 1 is perfectly cylindrical and V p is identical to V p . 
Since meniscus height variations yield radius variations (see equation (1)), achieving a 
relatively small fluctuation in crystal diameter is coupled to achieving a relatively small 
fluctuation in growth rate (see equation (2)). In the present embodiment, information on 
20 growth rate is stored in the first derivative of the meniscus height h(t) and, thus, in the 

second derivative of r(t). It is to be understood that a conventional PID controller ignores 
this basic quantity. 
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As described above, the ability to grow silicon crystals with controlled defect 
properties strongly depends on the ability to control growth rate. A key parameter in 
modern crystal growth is V /G, where G is the average axial temperature gradient at the 
melt/solid interface. Successful defect-free crystal growth is accomplished when the 
growth rate deviates from a critical value (e.g., VJG of about 0. 13 mnr/K min.) by less 
than 5-10%, according to hot zone design, for the entire length of crystal 31 . 
Advantageously, control circuit 51 controls the growth rate to a desired accuracy during 
the entire CZ crystal growth process even though the thermal conditions change 
continuously and specifically for particular hot zone configurations or process conditions. 
In this embodiment, control circuit 51 maintains diameter within a relatively low tolerance 
from the set value (e.g., 0.5 mm) and the growth rate (averaged over the control cycle) 
within a relatively low tolerance from the set value (e.g., 0.02 mm/min.) to grow perfect 
silicon. 

As shown by equations (1 ) and (2), variations in the diameter are a manifestation of 
the meniscus height oscillations around the steady-state value h s . Accordingly, it implies 
that the pull rate V p is not identical to the growth rate V u . Therefore, the requirement of 
small diameter deviations is actually coupled to the requirement of small growth rate 
deviations relative to the steady-state value, V . The growth rate V gs for the steady-state 
meniscus height h s is of principal importance for grow th control. Generally, V g (for any 
meniscus height) is expressed from the heat flux balance at the solid/melt interface by 
equation (3), below: 

V,(t) = [k s G s (t) - k m G m (t)] / Q (3) 

where G s (t) is the solid state temperature gradient (°C/mm); G ni (t) is the liquid state 
temperature gradient (°C/mm); k s is the solid phase Si conductivity (J/mm min K); k m is 
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the liquid phase Si conductivity (J/mm min K); and Q is the heat of fusion (J/mm 3 ). In this 
instance, k s = 1 .5 J/mm min K, k m = 3 J/mm min K and Q = 4. 1 J/mm 3 . 

The melt gradient G m is also sensitive to the meniscus height. In particular, G m is 
roughly inversely proportional to h. According to the invention, the derivative of V tt with 
5 respect to h provides an important control parameter, namely, the height response 

coefficient A h defined by equation (4), below: 

A h = dVJdh > 0 (4) 

In equation (4), the sign is positive since an increased height means a decreased gradient 
G n , 

10 If puller 13 pulls crystal 31 at V p = V gs , then the meniscus height will automatically 

go to the value of h s , according to equation (2), after a transient time of 1/A h (e.g., about 5 
minutes). After that, the diameter of crystal 3 1 automatically becomes constant while the 
actual growth rate V„ becomes identical to V . According to the invention, the process 
control executed by control unit 51 deduces from an observed diameter evolution and 

1 5 controls the pull rate at V p = V„ s . Next, a proper power increment brings the value of V gs 

to the set value of V sct . Increasing the power leads to a higher melt temperature below 
crystal 31 and, accordingly, to a higher G m (and to lower V ). The derivative of V„ s with 
respect to the power P is another important control parameter defined by equation (5), 
below: 



20 A p = dVJdP < 0 (5) 



Since thermal conditions change continuously during CZ crystal growth, there is 
the need to accommodate for changes in during the run, by intervening on the power 
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that power supply 27 provides to heater 21 so as to keep V gs close to the set growth rate 
V sct (t). The required power increment is defined by equation (6), below: 

6P = ~(V ;JS - V s ,)/A p (6) 

Applying such a power increment causes V to move toward the set value of V sc , after a 
5 transient time x corresponding to the 'crucible heat inertia'. For example, z is on the order 

of 20 minutes. 

Referring now to FIG. 3, the present invention controls the pull rate and heater 
power to provide growth rate control that is directly translated into simultaneous control of 
both diameter and growth rate. FIG. 3 illustrates the operation of control unit 51 according 

10 to a preferred embodiment of the invention for accomplishing such growth rate control. 

Beginning at step 69, control unit 51 records the radius function r(t) for a predetermined 
interval or period of time. Control unit 51 then processes the information at step 71. In 
processing r(t), control unit 51 first deduces current values of r, h and V at step 73. 
Appendix A describes a best fit analysis of the observed r(t) function for deducing these 

1 5 values according to the invention. 

In a preferred embodiment of the invention, the model-based controller embodied 
by control unit 51 processes the open loop radius variation r(t) recorded for the given 
interval of time to deduce the current values of r, h and V (and then V ). Proceeding to 
step 75, control unit 51 further processes the deduced values of r, h and V g to deduce a 

20 current value of the steady-state growth rate V us . At step 77, control unit 51 generates an 

output signal representative of an incremental change in pull rate as a function of the 
deduced variable V es . Control unit 51 then applies each pull rate increment, or pull action, 
to crystal drive unit 47 with a given periodicity (i.e., at the end of each observation 
interval) at step 79. FIG. 4 is an exemplary graph illustrating deviations in the radius of 
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crystal 31 relative to time and comparing the actual radius that results from an incremental 
change in pull rate to the desired radius at the end of the particular interval. 

Referring again to FIG. 3, control unit 51 observes r(t) over N observation periods, 
each of approximately equal duration. The period duration t p is preferably on the order of 
5 the meniscus response time 1/A h or somewhat longer (e.g., about 10 to 20 minutes). At 

step 81, after N periods and, thus, after N pull actions, control unit 51 accumulates N 
values of V gs . The current value of V gs (used in the power action recipe equation (6)) is 
deduced using the N accumulated values for V gs . Control unit 51 accumulates the values 
of V us at step 81 and deduces the current value of V gs for the power action accordingly at 

10 step 83. This accumulation is important to reduce an error in V gs caused by a noise in the 

measured diameter. At step 85, control unit 51 generates an output signal representative of 
an incremental change in heater power as a function of the deduced variables and then, at 
step 87, applies a power action to heater power supply 27 according to the power 
increment specified in equation (6). 

1 5 In one embodiment, the number of pull periods N that precedes the power action is 

about six so that the period between power actions is about one hour or more. Crystal 3 1 
grows about 2 cm in length during each period (for growth conditions required for perfect 
silicon). The output increments in pull rate and heater power are intended to minimize, 
respectively, diameter and growth rate deviations. 

20 The approach is thus to leave crystal 31 c at peace' for some time to obtain 

sufficient information on the r(t) function at an essentially constant pull rate. Then, the 
current variables r and h and the current values of V g and V gs can be reliably deduced even 
in the presence of measurement noise. These variables being know n, one can predict the 
future behavior of both r(t) and h(t) for any new value of V and, thus, choose the best new 

25 V p to reduce both the radius deviation and the height deviation. The second goal is even 

more important since the derivative dr/dt is then reduced (see equation ( 1 )). Therefore, a 
'calm behavior 1 of diameter during future periods is ensured. 
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FIG. 5 is an exemplary graph illustrating crystal radius deviations as affected by 
incremental changes in pull rate and heater power as controlled by the apparatus of FIG. I . 
FIG. 5(a) shows the difference between the actual radius r and a desired radius r sct over 
time. FIGS. 5(b) and 5(c) show the incremental pull and power actions, respectively. In 
5 the example of FIG. 5(c), control unit 51 applies the first power action after the first 

observation period of 10 minutes and, afterwards, applies such actions after every N 
periods where N = 6. 

Appendix B describes an exemplary pull rate recipe according to the invention. 

In view of the above, it will be seen that the several objects of the invention are 
10 achieved and other advantageous results attained. 

As various changes could be made in the above constructions and methods without 
departing from the scope of the invention, it is intended that all matter contained in the 
above description or shown in the accompanying drawings shall be interpreted as 
illustrative and not in a limiting sense. 
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APPENDIX A 

According to the invention, an analysis of the observed r(t) function yields the 
current values of r, h and V,, (and then VJ. For this purpose, the growth rate V g (t) in the 
basic equations (1) and (2) is expanded with respect to the meniscus height: 

5 V 8 = V gr +A h (h- h r ) (A-l) 

where h t and V, t are the current values (corresponding to the end of the observation 
period). As such, equation ( 1 ) is linear but equation (2) is essentially non-linear. Yet the 
solution for r(t) at specified values of h t , V uf and r f (the radius at the end of the observation) 
is easily obtained in the analytical form: 

1 0 r(t) = r f + A[ Y(h f - h s )x + (Y + h f - h s )Z(x + 1 - e x ) + Z 2 (x + 1 .5 - 2e x + 0.5e 2x )] (A-2) 

where x = A h t; and the time t is counted back from the end to the start of the observation 
period. 

The parameters Y and Z are growth/pull rate combinations: 

Y = V gr /A h 

15 Z-(V p - V gl )/A h 

and the capillarity parameter is defined by: 



A = 2h s / a 2 
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The theoretical function r(t) includes the current system variables r t and h t , the 
current growth rate V ut - and the basic system parameter A h (i.e., the height response 
coefficient), all of which are to be determined. In a preferred embodiment, these four 
numbers are fitted to the measured function r exp (t) to get the least square deviation of r(t) 
5 from r (t). It is expected that noise in the function r cxp (t) prevents it from providing a 

reliable value for A h . For this reason, the present invention specifies a value for A h , which 
in turn permits control unit 5 1 to reliably deduce the three numbers r i? h, and V f by the 
best fit routine described above. Preferably, the present invention also recalculates V„ s , the 
growth rate at the steady-state meniscus height, which is a relevant quantity defining a 
10 proper intervention on the pull rate and power. This important quantity is recalculated 

from the fitted numbers for V Kt - and h t : 



V, s = V^ f - A h (h r -h,) 



(A-3) 
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APPENDIX B 



At the end of the observation period, the current values of r, h and V us are 
determined, obtained by the best fit routine described above in Appendix A. The values r t 
and h, are the final ones for the previous observation period but are starting values for the 
5 next period. Accordingly, these values are now denoted r G and h G . Integration of the 

dynamics equations (1 ) and (2) is now made in the forward direction, from the start to the 
end, using the expansion of V (J (h) around the steady-state height. In this manner, the 
predicted new final values for h and r are obtained: 

h = h s + Z + (h 0 - h s )e * (B-l) 

10 r = r tl - A[YZx 4- (2Z + Y)(h 0 - h s - Z)(l - e x ) + 0.5 (h u - h s ) 2 (l - e" 2x )] (B-2) 

where t p is the period duration; and x = A h t p . The parameters A, Y and Z are as defined 
above except V u . replaces V uf . 

To select the new value of the pull rate (entering the equations through the 
parameter Z = (V - V )/A h ), the current values of r G -r sct and h tl -h s are first checked. 
15 Preferably, the values are 'safe' to the prescribed diameter deviation (e.g., the radius 

deviation should not exceed some limit, Ar = 0.5 mm). The 'safe' situation means that 
pulling at the rate of V = will not cause the radius to deviate by more than Ar. From 
equation (B-2), the predicted r for this case, where Z = 0, is: 



r = r„ - A[ Y(h„ - h s )( 1 - e *) + 0.5 (h 0 h s ) 2 ( 1 - c 2 ^)] 



(B-3) 
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So if |r-r ( J is within the Ar limit, then the recipe for the pull rate is simplified: 

V p = V gs (B-4) 

As noted above, pulling at this rate automatically brings the meniscus height to the desired 
steady-state value of h s since h-h s =(h 0 -h s )e" x according to equation (B-l ) However, 
noise-induced error in the deduced value of V„ s , for example, may cause the predicted 
radius from equation (B-3) to go beyond the prescribed error of Ar. In such an instance, a 
proper height deviation may be introduced to reduce the radius deviation. For this 
purpose, V p should deviate from V us . The proper V p should be such that the new final 
values of r-r sel and h-h s are of the same sign (e.g., positive deviation r-r sct will decay at 
positive h h s according to equation (2)). The fastest subsequent decay is achieved if: 

(r-rJ/(h-h l ) = AV p /A h (B-5) 

according to equations (B-l) and (B-2). 

Therefore, V p is chosen to bring r and h to the desired condition as specified in 
equation (B-5). For this purpose, a computational cycle over V p is made to look at the 
final r and h values specified as the functions of V p by equations (B-l) and (B-2) defining 
h and r, respectively. The value of V p resulting in the relation of equation (B-5) is selected 
as the recipe for the pull action. In one preferred embodiment, limitations are imposed on 
V p (e.g., to be positive but less than a threshold) and on |h-h s ! (e.g., to be less than or 
equal to about 1 mm) to ensure the applicability of the basic meniscus equation (1). Thus, 
the finally accepted V is the value that brings the ratio (r r sct )/(h-h s ) as close as possible 
to the goal ratio AV /A h of equation (B-5) w ithin the imposed limitations. 



